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ABSTRACT 



Aims. We study the Be star S Cen circumstellar disk using long-baseline interferometry which is the only observing technique capable 

of resolving spatially and spectroscopically objects smaller than 5 mas in the H and K band. 

Methods. We used the VLTI/AMBER instrument on January 6, 8, and 9, 2008, in the H and K bands to complete low (35) and medium 

(1500) spectral resolution observations. 

Results. We detected an oscillation in the visibility curve plotted as a function of the spatial frequency which is a clear signature of a 

companion around S Cen. Our best-fit soltution infers a binary separation of 68.7 mas, a companion flux contribution in the K band of 

about 7% of the total flux, a PA of 1 17.5 °, and an envelope flux around the Be primary that contributes up to about 50 % of the total 

flux, in agreement with our Spectral Energy Distribution (SED) fit. The envelope size is estimated to be 1.6 mas in K but no departure 

from spherical symmetry is detected. 
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1. Introduction 

We have initiated an observational campaign to determine the 

global parameters of Be stars and the origin of the so-called "Be- 

phenomenon" using the VLTI interferometer, which combines, 

for the first time, milli-arcsecond spatial resolution with good 

spectral resolution of up to 12000 in the near-IR Kband. We are 

therefore able to study the geometry and kinematics within the 

, ^ circumstellar environment. Three Be stars have been observed 

K^ , with the VLTI instruments and the global physical properties of 

5_i ■ these stars have been found to vary strongly from one star to 

Cd ' another. 

Achernar was the first Be star observed with the VLTI us- 
ing the VINCI beam recombiner. These observations exhib- 
ited a strongly flattened star due to an almost critical rotation 
(Domiciano de Souza & Kervella 2003) and an extended polar 
wind that contributed to almost 4% of the total flux in the H 
band (Kervella & Domiciano de Souza 2007) with no detectable 
equatorial disk. Nevertheless, its circumstellar environment was 
known to be highly variable with a quasi-cyclic formation and 
a dissipation timescale of about 13 years proposed by Vinicius 
et al (2007). Kanaan et al. (2008) modeled these variations suc- 
cessfully as a possible outburst of circumstellar matter that has 
been ejected and was propagating within the stellar environment 
with an expansion velocity of 0.2km s . 

a Arae was the second star observed with the VLTI, using 
both MIDI and AMBER instruments, and the observational re- 
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suits were published by Chesneau et al. (2005) and Meilland 
et al. (2007 a). The data interpretation was completed with the 
SIMECA code developed by Stee et al. (1994), which is able 
to simulate the gaseous circumstellar environment of active hot 
stars. Two distinct regions in the circumstellar environment of 
a Arae are clearly present : a thin Keplerian rotating equato- 
rial disk and a rapidly expanding enhanced polar wind. As for 
Achernar, the stellar rotational velocity, which was readily mea- 
surable due to the flatness of the projected envelope onto the sky 
plane, was found to be almost critical. By combining MIDI and 
AMBER data, they found that the disk size did not vary as a 
function of the observing wavelength, which was a strong indi- 
cation of possible disk truncation, as confirmed by more "classi- 
cal" spectroscopic measurements. 

Finally, the third Be star studied using the VLTI instruments, 
was k CMa, observed with AMBER in MR mode and presented 
in Meilland et al. (2007b). For this star, a non-Keplerian rotat- 
ing disk was detected. The star was found to rotate far below 
its critical velocity, which appears to indicate that critical rota- 
tion may not be a requirement of the "Be phenomenon". Since 
it has an early spectral type of B2 V, the radiation pressure might 
however be sufficient to initiate a mass loss and possible disk for- 
mation. This hypothesis agrees with the Cranmer (2005) statis- 
tical study. Finally, the spectro-differential interferometric mea- 
surements were able to detect the presence of an inhomogeneity 
in the circumstellar envelope of this star, which could be inter- 
preted in terms of the "one-armed" viscous disk framework of 
Okazaki( [T997j l. 

To increase the number of Be stars observed by the VLTI and 
initiale a statistical study of the geometry and kinematics of their 
surroundings, and study the possible dependence of the "Be phe- 
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nomenon" on the spectral type of the central star, we requested 
new AMBER observations of 6 close and bright enough Be stars 
of spectral types ranging from Bl to B8 of which 5 Cen was one 
star. 

This letter is organized as follows: in Sect. 2 we briefly intro- 
duce the Be star 5 Cen and its fundamental parameters, in Sect. 
3 we present our VLTI/ AMBER observations and the data re- 
duction, and in the last section our results are summarized and 
conclusions are presented. 
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Fig. 1. Left : Names, lengths, and RA. for each baseline of the 
two triplets of observation. Note that the two measurements from 
January 8 with A0-K0-G1 are merged into Bi, B2, and B3 base- 
lines, since the projected baseline length and RA. does not vary 
significantly between them. B4, B 5 , and B 6 baselines correspond 
to the E0-G0-H0 triplets. Right: (u,v) plane coverage obtained 
during the 5 Cen observing run. Fortunately, the largest projected 
baseline is aligned with the shorter ones. 

2. 5 Cen 

5 Cen (HR 4621, HD 105435, HIP 59196) is usually classified as 
a B2IVne star. Levenhagen & Leister (2006) considered a spec- 
troscopic analysis that did not take account of the effects of grav- 
ity darkening and proposed a B2Vne classification. It is a vari- 
able star with a brightness variation between magnitude 2.65 and 
2.5 1 . 6 Cen has the same proper motion as the nearby stars HD 
105382 (which is also a Be star) and HD 105383 and may there- 
fore be member of a triple visual star system. A precise radial 
velocity (RV) measurement campaign of B stars in the Scorpius- 
Centaurus association conducted by Jilinski et al. (2006) mea- 
sured a RV for 5 Cen of 3.8+ 2.8 km s . The RV determined by 
various authors from the SIMBAD database also exhibits a large 
dispersion, which might indicate possible binarity, as we see in 
the following. On the other side, Mason et al. (1996) and Mason 
(2008), using U.S. Naval Observatory (USNO) and Center for 
High Resolution Astronomy (CHARA) speckles cameras did not 
detect any binarity for this star. 

Based on an autocorrelation analysis of Hipparcos photom- 
etry, Percy et al. (2002) presented a self-correlation diagram 
with strong minima at 0.55 and 1.1 days, in contrast, Balona 
(1990) reported a tentative period of 1.923 days, implying that 

5 Cen is a short-period variable Be star. Studying line profile 
variations in 1999 HEROS spectroscopic data, Rivinius et al. 
(2003), confirmed both values at 0.532d and 1.1 39d but with 
a formally somewhat higher significance of the longer period 
which was adopted in their paper. They concluded that 5 Cen 
line profile variability was modeled successfully as a non-radial 
pulsator with I - m - +2, as for the Be star to CMa, but with a 
strongly asymmetric power distribution of the line profile vari- 
ability with a unidentified cause. The star's distance was esti- 
mated by Hipparcos to be 121 + 1 lpc, and its Kand V magnitude 
to be 2.7 and 2.56, respectively . Lenorzer et al. (2002) identified 

6 Cen as a B star with emission lines observed in the ISO Short 
Wavelength Spectrometer Post-Helium program of the Infrared 



Obs. time 


Target 


Triplet 


Mode 


Seeing 


06/01/08 08:59 


6 Cen 


A0-K0-G1 


MR 


1.02" 


06/01/08 09:11 


HD109538 


A0-K0-G1 


MR 


1.13" 


06/01/08 09:26 


6 Cen 


A0-K0-G1 


MR 


1.19" 


08/01/08 08:48 


HD90798 


E0-G0-H0 


LR 


0.50" 


08/01/08 09:11 


6 Cen 


E0-G0-H0 


LR 


0.94" 


09/01/08 08:07 


HD47536 


E0-G0-H0 


MR 


0.68" 


09/01/08 08:31 


6 Cen 


E0-G0-H0 


MR 


0.74" 


09/01/08 08:50 


HD109538 


E0-G0-H0 


MR 


0.87" 


09/01/08 09:08 


S Cen 


E0-G0-H0 


MR 


1.22" 


09/01/08 09:28 


HD109538 


E0-G0-H0 


MR 


0.71" 



Table 1. Observations log of 5 Cen and its calibrators during the 
January 2008 Be AMBER GTO run. 



Space observatory (ISO). The fundamental parameters deter- 
mined for rapidly rotating B-type stars by Fremat et al. (2005), 
taking into account gravitational darkening were T e //=22384 
+ 446 K, log g=3.942 ±0.053, and Vsim=263±14 km s _1 . 
Its linear critical equatorial rotation velocity was V e =527+29 
km s and its inclination angle estimated to be z'=41.6+2.3°. 
Levenhagen & Leister (2004) estimate its age to be log Age = 
7.20 +0. 15 yr, its luminosity log L/Lq=3.70 +0. 10, and its mass 
8.6 ±0.3 M Q . 



3. Observations and data reduction 

5 Cen was observed with AMBER in low and medium spec- 
tral resolution modes (spectral resolving powers of R=35 and 
R=1500, respectively) during the observing runs 080.C-0099 
and 080.D-0131, completed in January, 6, 8, and 9, 2008. 
Five measurments were obtained with the Auxilliary Telescopes 
(1.8m diameter) on the A0-K0-G1 and E0-G0-H0 baseline 
triplets. HD 109593 (apparent diameter of 1 .638± 0.021mas in 
the CHARM2 catalog from Richichi et al. [2005b . HD 47536 
(apparent diameter of 1.69±0.03mas, also in CHARM2), and 
HD 90798 (apparent diameter 1.70±0.51 mas using the getCafl 
tool) were also observed to calibrate 5 Cen visibilities. The log 
of these observations are presented Table 1 . 

The data were reduced using the AMBER data reduction 
software developpedby the JMMC: amdlib, version 2. 1 . It com- 
putes a fringe fitting of the AMBER data, after correcting for 
detector bias and flat-field (Millour et al., 2004 and Tatulli et al., 
2007). Following that step, we adopted two different approaches 
to calibrate the data: 

- The "standard" approach, which retains 20% of the total 
number of frames according to the fringe SNR. This allows 
us to limit the fringe smearing effect (also called "jitter" ef- 
fect) and the coherence length visibility loss due to the at- 
mospheric OPD (described in Millour et al. 2007). The fol- 
lowing steps (transfer function estimation, calibration) were 
completed using a new set of custom-made scripts described 
by Millour et al. (f20083). 

- An approach using bootstrapping, roughly explained in 
Tatulli et al., (2007) and described in more detail in 
Cruzalebes et al. d2007l l. A frame-selection scheme that re- 
jects frames with a negative flux, a high OPD variation from 
frame to frame, and a SNR of less than 1, is also applied. The 
calibration procedure is then completed the same way as for 
the previous method. 



1 available at http : //mscweb . ipac . caltech.edu/gcWeb 

2 The scripts are available to the community at the following web- 
page: http://www.mpifr-bonn.mpg.de/staff/fmillour 
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5 Cen SED 




Fig. 2. 6 Cen Spectral Energy Distribution (SED) reconstructed 
from various sources in the literature : Morel (1978), 2MASS, 
IRAS magnitude and flux measurements (diamonds), IUE spec- 
tra (solid line in the UV), and ISO spectra (solid line in the 
IR). The dashed line is a Kurucz model for T e ff - 21000K, 
R* = 5.9R Q , log g = 4, and d = 121pc. 

The two methods provide the same results within the error 
bars and we therefore choose to use the calibrated data from 
the first method later in the paper. Since the aim of the paper 
is to determine the binary parameters, we exclude the January 
9 MR data for which projected baselines are similar to the 
January 8 LR ones and the data quality is significantly lower. 
Finally, we merged the two January 5 MR measurements ob- 
tained with almost the same projected baseline, and determined 
the wavelength-dependent visibilities for the six baselines pre- 
sented Fig. 1. The corresponding (u,v) plane coverage of the sci- 
ence target is also plotted in Fig 1 . 

4. Results 

For a single Be star, assuming that the visibility we measure, V, 
is only due to the central star and its circumstellar disk we can 
write: 



V = 



V em ,F em , + V+F+ 



(1) 



where V em , and F em , are respectively the envelope visibility and 
flux in the continuum , V* and F* the star visibility and flux in 
the continuum and F to , = F cnv +F*. To measure accurately the ex- 
tension of a circumstellar envelope from an interferometric mea- 
surement, we therefore need to determine both F* and F em , at the 
observing wavelength. This can be achieved directly by analyz- 
ing data acquired by interferometry if, at least, two baselines 
resolve the envelope completely but not the star (i.e, V fm ,=0 and 
V* ~ 1 ). In the case of our 5 Cen AMBER data, even the longest 
baseline (i.e, B3 = 127.9m) is too short to fully resolve the cir- 
cumstellar envelope. 

To determine F* and ¥ env in H and K bands to the highest 
possible accuracy, we started our modeling by constructing the 
6 Cen Spectal Energy Distribution (SED) using various sources 
in the literature (see Fig. 2). A Kurucz model, fitted from UV and 
visible measurements, for which the circumstellar contribution 
was negligible, was used to determine the stellar contribution to 
the flux for wavelengths ranging from UV to Far-IR. The values 
of the best-fit Kurucz model parameters, which were consistent 
with those found in the literature and were introduced in Sect. 
1, are presented Table 3. This method allows us to characterize 
precisely the circumstellar IR-excess. Taking into account the 
uncertainty in the spectral class of 6 Cen and its small variability 



in the near IR, we obtain F em ,= 0.52+0.06 F„„ in the K band and 
0.38+0.08 F lol in the H band. 

K and H band calibrated visibilities obtained with the small 
aligned triplet in LR mode (i.e. B4, B5, and Be) are plotted as 
a function of the spatial frequency in Fig. 3. These data exhibit 
clearly a sinusoidal modulation with a period of T4,5,6=22 cy- 
cles/arcsec and an amplitude of A456=0.15±0.02. This varia- 
tion cannot be reproduced by a simple star + circumstellar en- 
velope model regardless of the geometry of this second com- 
ponent (e.g. uniform or Gaussian disk, or ring), and is usually 
produced by a companion. The separation between the two com- 
ponents in the direction of the projected baselines (i.e., -11 1°) is 
given by 84,5,6= l/To=44.8±0. 5 mas, and its fluxbyF f =A4,5,6/2 = 
0.075 ±0.0 lF fo ,. We note that this value may be biased towards a 
dimmer flux ratio since the separation is ~25% of the telescopes 
PSF (~ 240mas with the ATs at 2.1/an) and the AMBER instru- 
ment uses monomode fibers (Tatulli et al. 2004). However, it is 
beyond the scope of this paper to correct for this effect. 

To reproduce the small slope in the visibilty as a function 
of spatial frequency, we need to introduce a small circumstel- 
lar envelope (i.e., D=2±0.5mas) around one of the two compo- 
nents. We note that the data also indicate some small differences 
in modulation amplitude between the H and K band data. The 
characterization of these differences, probably due to the varia- 
tion in the binary flux ratio between these spectral bands, is again 
beyond the scope of this letter and is not described here. 

The Bi, B2, and B3 data are also plotted in Fig 3. These 
measurements place additional constraints on the binary and the 
circumstellar envelope. The B2 and B3 data show a sinusoidal 
modulation with periods corresponding to a separation between 
the binary components in the projected baseline directions of 
S2=62+5 mas and S3=52±6 mas, respectively. The amplitude of 
the modulation is A2=A3=0. 15+0.05 which is in agreement with 
the measurement for the January 8 LR data (i.e., Ao=0. 15+0.02). 
The separation measured by B3 is consistent with that measured 
for data acquired with the LR triplet; this is a good indication that 
our measurements are not biased by some instrumental effects 
since B3 is almost aligned with the LR triplet (i.e. P.A.=-1 16° for 
B3, and P.A.=-1 1 1° for B4, B5, and Be). The difference between 
the two measurements can be due either to a inaccurate spec- 
tral calibration of AMBER LR data (up to 0. l/im) or to a small 
change in the projected separation due to the 5" of variation in 
RA. between the two data sets. Finally, there is no evidence of a 
sinusoidal modulation in the Bi data. Since the modulation may 
not be detected if its period is at least 3 times the observation 
range in spatial frequency (i.e., T>50 cycles/arcsec), the binary 
separation along this projected angle, Si, is lower than 20 mas. 
The average visibility of each baseline (i.e. Vi=0.82, ^2=0. 86, 
and V3=0.8 1) indicates clearly the presence of an extended com- 
ponent, probably a circumstellar envelope with an extension of 
about 2±0.5mas. 

By measuring the projected separation in at least two dif- 
ferent directions in the sky, we can determine the modulus and 
the orientation of the binary separation on the sky plane. Our 
most appropriate binary star + circumstellar envelope model is 
obtained for a set of values given in Table 2. The modeled visi- 
bility obtained in K band for each of the six baselines is plotted 
with the observed measurements in Fig. 3. We note that to ver- 
ify that the binary signal did not originate in the calibrators we 
checked that the modulation was evident in the science object 
raw visibilities and not on the calibrator ones. 

Since AMBER operates simultaneously with three tele- 
scopes, we are able to measure the closure phase for each triplet. 
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Fig. 3. Left: LR K (dark) and H (gray) bands visibilities for the January, 8 2008 baselines (B4, B5, and B^) as a function of the spatial 
frequency. Second picture from the left to the right: MR visibilities for the January, 5 2008 baselines (Bi, B2, and B3) as a function 
of the spatial frequency. For all figures, measurements are given by the error bars, whereas the solid line corresponds to our best-fit 
binary + uniform disk circumstellar envelope modeled using parameters from Table 2. 



Such measurements place additional constraints on our models. 
Figure 4 shows a reasonable agreement between observed and 
modeled closure phase using parameters from Table 2. 

Finally, using R*=5.9R Q and d=121pc, the 68.7mas sep- 
aration between 6 Cen and its companion corresponds to 
300+50R*. Assuming a circular orbit, this provides a lower limit 
to the semi-major axis. The companion spectral class can be 
inferred using both 6 Cen spectral class and K band flux ratio 
determined from the amplitude of the modulation in the inter- 
ferometric data. The companion spectral class should range be- 
tween B4V and A0III, and its mass between 4 and 7 M Q . Using 
Kepler's third law of planetary motion, P = -\Ja 3 /M tot , where 
P is the period in years, a is the semi-major axis in AU, and 
M tot is the binary system total mass in M , we can determine 
a lower limit to its period, i.e. 4.6 years with M, ,=15.2M and 
a=250R*=6.9AU. 

Additional interferometric observations, including long-term 
monitoring, would be required to determine why these results 
disagree with those of the speckles observations of Mason (1997, 
2008) and to fully determine the projected orbit of the compan- 
ion of 6 Cen. Spectroscopic follow-up should also be completed 
even though the Doppler shifts due to binarity are difficult to 
measure if the star is a rapid rotator (i.e, V sin i= 263km s -1 ), 
and the profiles are affected by non-radial pulsations. The cir- 
cumstellar environment geometry and kinematics will be stud- 
ied in a forthcoming paper, which will include all January 2008 
VLTI/ AMBER medium resolution differential visibilties as well 
as VLTI/MIDI data. 



Distance 


121 ±llpc 


T'cff 


21000 ± 1000 K 


Stellar Radius 


5.9 ±0.5 R© 


l°g£ 


3.75 ± 0.25 


Companion flux (K band) 


0.07±0.0.1 F fof 


Separation 


68.7±0.5 mas 


Position Angle (PA) 


117.5±0.5° 


Envelope flux (K band) 


0.52+0.06 F M 


Envelope diameter (K band) 


1.6±0.4mas 



Table 2. Parameters and results for the binary system + uniform 
disk circumstellar envelope obtained for 5 Cen. Note that there is 
a 180° ambiguity in the position angle since the AMBER closure 
phase data reduction process is not fully validated. 
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Fig. 4. Observed and modeled closure phase for the B4B5B 6 
(left) and B1B2B3 (right) triplets plotted as a function of the av- 
erage spatial frequency. For all figures, the measurements are 
represented by the error bars whereas the solid line corresponds 
to our best binary + uniform disk circumstellar envelope mod- 
eled using parameters from Table 2. 
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